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Abstract 

Econophysics has developed as a research field that applies the for¬ 
malism of Statistical Mechanics and Quantum Mechanics to address 
Economics and Finance problems. The branch of Econophysics that 
applies of Quantum Theory to Economics and Finance is called Quan¬ 
tum Econophysics. In Finance, Quantum Econophysics’ contributions 
have ranged from option pricing to market dynamics modeling, behav¬ 
ioral finance and applications of Game Theory, integrating the empir¬ 
ical finding, from human decision analysis, that shows that nonlinear 
update rules in probabilities, leading to non-additive decision weights, 
can be computationally approached from quantum computation, with 
resulting quantum interference terms explaining the non-additive prob¬ 
abilities. The current work draws on these results to introduce new 
tools from Quantum Artificial Intelligence, namely Quantum Artifi¬ 
cial Neural Networks as a way to build and simulate financial market 
models with adaptive selection of trading rules, leading to turbulence 
and excess kurtosis in the returns distributions for a wide range of 
parameters. 

Keywords: Finance, Econophysics, Quantum Artificial Neural 
Networks, Quantum Stochastic Processes, Cognitive Science 
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1 Introduction 


One of the major problems of financial modeling has been to address com¬ 
plex hnancial returns dynamics, in particular, excess kurtosis and volatility- 
related turbulence which lead to statistically signihcant deviations from the 
Gaussian random walk model worked in traditional Financial Theory (Arthur 
et ah, 1997; Voit, 2001; Ilinsky, 2001; Focardi and Fabozzi, 2004). A main 
contribution of Econophysics to Finance has been to address this problem us¬ 
ing the tools from Statistical Mechanics and Quantum Mechanics, within the 
paradigmatic basis of Systems Science and Complexity Sciences (Anderson 
et ah, 1988; Arthur et ah, 1997; Voit, 2001; Ehrentreich, 2008). 

Econophysics is currently a major research area that has combined inter- 
disciplinarily Finance and Economics, Complex Systems Science, Statistical 
Mechanics, Quantum Mechanics and Cognitive Science to address notions 
and problems in Economics and Finance (Anderson et ah, 1988; Arthur et 
ah, 1997; Voit, 2001; Brunn, 2006; Ehrentreich, 2008; Piotrowski and Slad- 
kowski, 2001, 2002, 2008; Saptsin and Soloviev, 2009, 2011). 

There are two major branches in Econophysics: Classical Econophysics 
(based on Classical Mechanics) and Quantum Econophysics (based on Quan¬ 
tum Mechanics). In Finance, Quantum Econophysics has been applied to 
option pricing (Segal and Segal, 1998; Baaquie et ah, 2000; Baaquie and 
Marakani, 2001; Baaquie, 2004; Baaquie and Pan, 2011), hnancial turbulence 
modeling (Gongalves, 2011, 2013) and as an approach to the formulation of 
hnancial theory, regarding price formation and basic market relations (Pi¬ 
otrowski and Sladkowski, 2001, 2002, 2008; Khrennikov, 2010; Haven and 
Khrennikov, 2013; Gongalves, 2011, 2013). Choustova (2007a,b), in partic¬ 
ular, argued for the introduction of a quantum-based approach to Financial 
Theory as a way to incorporate market cognition dynamics in hnancial price 
formation. 

The quantum-based approach goes, however, beyond a good match to 
price dynamics and turbulence modeling. The growing empirical evidence of 
quantum interference signatures in human cognition, when faced with deci¬ 
sion problems, has led to the development of a Quantum Theory-based Cog¬ 
nitive Science forming a theoretical ground for Econophysics modeling, with 
strong implications for Finance (Busemeyer and Franco, 2010; Busemeyer 
and Bruza, 2012; Wang and Busemeyer, 2013; Busemeyer and Wang, 2014; 
Khrennikov, 2010; Haven and Khrennikov, 2013; Zuo, 2014; Khrennikov and 
Basieva, 2014). 
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The main research problem regarding Quantum Theory-based Cognitive 
Science applied to Finance can be expressed as follows; if there is emprical 
support to the fact that human cognition, in decision problems, leads to a 
decision behavior computationally isomorphic to quantum adaptive compu¬ 
tation (Busemeyer and Franco, 2010; Busemeyer and Bruza, 2012; Wang and 
Busemeyer, 2013; Busemeyer and Wang, 2014; Zuo, 2014; Khrennikov and 
Basieva, 2014; Gongalves, 2015), then, the modeling of hnancial market dy¬ 
namics needs to work with models of behavior that incorporate, in their prob¬ 
abilistic description, quantum interference terms (Khrennikov, 2010; Haven 
and Khrennikov, 2013). 

This main research problem has led to the growth and development of 
research lines on Cognitive Science, working from Quantum Computer Sci¬ 
ence and Quantum Information Theory, with direct implications in Finance 
and Economics, supporting the expansion of Quantum Econophysics (Khren¬ 
nikov, 2010; Haven and Khrennikov, 2013), in particular, in regards to Fi¬ 
nance: opening up the way for research on Quantum Artihcial Intelligence 
(QuAI) applications to hnancial market modeling (Gongalves, 2011, 2013). 

The current work contributes to such research by introducing Quantum 
Artihcial Neural Networks (QuANNs) for hnancial market dynamics and 
volatility risk modeling. In particular, recurrent QuANNs are used to build 
a model of hnancial market dynamics that incorporates quantum interfer¬ 
ence and quantum adaptive computation in the probabilistic description of 
hnancial returns. The resulting model shows a quantum-based selection of 
adaptive rules with consequences for the market dynamics, leading to excess 
kurtosis and turbulence with clustering volatility, price jumps and statisti¬ 
cally signihcant deviations from Gaussian distributions, for a wide range of 
parameters. 

The work is divided in two parts that are developed in sections 2 and 
3. In section 2, the QuANN model is built, simulated and studied, while, 
in section 3, a rehection is provided on the possible role and contributions 
of QuAI applied to hnancial modeling. Regarding the main work, which 
is developed in section 2, the structure of this section is divided in three 
subsections. 

In subsection 2.1, we review a general framework for Glassical Econo¬ 
physics modeling of hnancial market price formation in which Farmer’s mar¬ 
ket making model (Farmer, 2002) is reviewed and combined with multiplica¬ 
tive components, namely: multiplicative volatility components and a mar¬ 
ket polarization component are introduced in the market making model and 
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linked to trading volume and bullish versus bearish polarization. 

In subsection 2.2, we introduce the general formalism of QuANNs, includ¬ 
ing main notions that form the groundwork for the hnancial market model. In 
subsection 2.3, we build the hnancial market model using a Quantum Neural 
Automaton (QuNA) structure and simulate the resulting artihcial hnancial 
market, addressing its main results in regards to turbulence and volatility 
risk, leading to statistically signihcant deviations from the Gaussian returns 
distribution. 

In section 3, the problem of deviations from the Gaussian random walk 
is addressed in its relation to Econophysics and Nonlinear Stochastic Models 
of market dynamics, allowing for a rehection on the possible contributions 
of QuAI and QuANNs for establishing a bridge between the evidence of 
quantum interference patterns observed in human decision making and a 
computational basis for nonlinear probability dynamics in Finance coming 
from a linear unitary evolution of networked quantum computation. 

2 A QuANN-Based Financial Market Model 

2.1 Price Formation and Financial Returns 

Following Farmer (2002) and Ilinski (2001), hnancial market price formation 
can be linked to unbalanced market orders M, where M > 0 corresponds to 
an excess demand while M < 0 to an excess supply, such that, for a hnancial 
risky asset, traded at discrete trading rounds of duration At, the asset price 
at t, S{t) depends upon the previous price S{t — At) and the market orders 
that arrive during the trading round. A few basic assumptions, in Glassical 
Econophysics, determine the structure for the relation between market orders 
and the new price (Farmer, 2002; Ilinsky, 2001); 

• The price is assumed as a hnite increasing function of the previous price 
and order size M{t): 


S{t) = fs{S{t-At),M{t)) (1) 

• If the order size is null M{t) = 0 the market clears for equal supply and 
demand, so that there is no market impact (the price stays unchanged): 

fsiSit-At),0) = Sit-At) ( 2 ) 
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• There are no arbitrage opportunities associated with a sequence of 
trades that sum zero (a repeated trading through a circuit); 

• Gauge invariance with respect to currency units, so that the only pos¬ 
sible combination for prices to enter is S{t)/S{t — At), such that; 


w^) ^ — w^) — ^ ^ 


(3) 


The result of these four assumptions is the general form for F in Eq. (3) given 
by (Farmer, 2002; Ilinsky, 2001); 

M(t) 

F{M{t)) = e— (4) 

where A is a liquidity parameter, also called market depth (Farmer, 2002). 
The result from Eq.(4), replaced in Eq.(3) is the following dynamical rule; 

S{t) = S{t — At)e^^ (5) 

or, taking the logarithms, the log-price rule (Farmer, 2002; Ilinski, 2001); 

In S{t) = In S{t — At)-\ -^ (6) 

There are two dynamical components to M{t): the sign, which can either 
be positive (excess of buy orders) or negative (excess of sell orders), and the 
volume of unbalanced market orders, which is linked to the order size. 

Within financial theory, the order size can be worked from a systemic 
market dynamics that leads to the formation of consensus clusters regarding 
the decision to invest greater or smaller amounts, or, alternatively, to sell 
greater or smaller amounts. The adaptive management of exposure to asset 
price fluctuation risk, on the part of market agents, given information that 
impacts asset value leads to a two-sided aspect of computation of hnancial 
information by the market system; on the one hand, there is the matter 
whether each new information is good [bullish) or bad [bearish), in terms of 
asset value, on the other hand, there is the degree to which new information 
supports the decision to buy or sell by different amounts (the market volume 
aspect). 

A social consensus dynamics coming from market computation can be 
linked to consensus clusters affecting the market unbalance, so that the pos¬ 
itive or negative sign can be addressed, within Econophysics, in terms of a 
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notion of spin. In Physics the spin is a fundamental degree of freedom of 
held quanta that behaves like angular momentum, the spin quantum num¬ 
bers assume integer and half-integer values, the most elementary case of half 
integer spin is the spin-1/2. 

Considering a three dimensional axes system, if a spin-1/2 particle’s spin 
state is measured along the ; 2 -axis then there are two fundamental orienta¬ 
tions spin up and spin down, in Complex Systems Science these two orienta¬ 
tions are assumed and worked mainly from the statistical mechanics of Ising 
systems as models of complex systems (Kauffman, 1993), which constituted 
early inspiration for Econophysics models of hnancial markets (Vaga, 1990; 
liori, 1999; Lux and Marches!, 1999; Voit, 2001). These models allowed for 
the study of polarization in market sentiment, working with the statistical 
mechanics of Ising systems, allowing direct connections to Cognitive Science 
(Voit, 2001). 

The market volume, on the other hand, has been addressed, within Fi¬ 
nancial Theory, by multiplicative processes (Mandelbrot, et ah, 1997; Man¬ 
delbrot, 1997), drawing upon Mandelbrot’s work on turbulence in Statistical 
Mechanics, as reviewed in Mandelbrot (1997). The multiplicative stochastic 
processes, worked by Mandelbrot and connected to multifractal geometry, led 
to Mandelbrot et al.’s (1997) Multifractal Model of Asset Returns (MMAR), 
which also inspired modihed versions using multiplicative stochastic processes 
with Markov switching in volatility components (Calvet and Fisher, 2004; 
Lux, 2008). 

Considering Eq.(6), a spin-1/2 like model can be integrated as a binary 
component in a multiplicative model that includes market volume, by way of 
a multiplictive decomposition of M{t) in a market polarization component 
a{t) = ±1, and N trading volume-related volatility components, so that we 
obtain; 


M(t )= (n j p) 

where each volatility component Vk{t) can assume one of two values Uq or 
Vi = 2 — vq. If 0 < Uq < 1; then vq corresponds to a low volatility state, 
while Vi to a high volatility state (vq diminishes the returns’ value while vi 
amplihes the returns value like a lever). The logarithmic returns for the risky 
asset, in this approach, are given by; 
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( 8 ) 


The binary structure assumed for the N components plus the market po¬ 
larization, makes this model a good starting point for QuANN applications, 
since QuANNs also work from a binary computational basis to address neural 
bring patterns. 

On the other hand, QuANNs open up the possibility for dealing with 
the multiplicative models in such a way that the probabilities, rather than 
being introduced from a top-down ex-ante bxed state-transition probability 
distribution, change from trading round to trading round, being the result of 
the quantum computational process introduced for each returns’ component. 

QuANNs also allow one to incorporate the empirical evidence that human 
cognition, when addressing decision between alternatives, follows a dynamics 
that is computationally isomorphic to quantum computation applied to De¬ 
cision Science, leading to interference effects with an expression in decision 
frequencies (probabilities), which means that, when considering probabili¬ 
ties for human behavior, the theoretical framework of networked quantum 
computation may be more appropriate for the dynamical modeling of human 
systems. 

In the quantum description, Eq.(8) will be expressed in operator form on 
an appropriate Hilbert space, with the returns operator eigenvalues being ad¬ 
dressed from the QuANN structure, which works with quantum hits {qubits), 
whose computational basis states decribe the neuron’s firing pattern in terms 
of firing (ON) and non-firing (OFF)Q In order to build the market model, 
however, we need to introduce, first, a general framework for QuANNs which 
will then be applied to the risky asset price dynamics modeling. 

2.2 Quantum Artificial Neural Networks 

The connection between Quantum Computer Science and ANNs has been 
object of research since the 1990s, in particular, in what regards quantum 
associative memory, quantum parallel processing, extension of classical ANN 
schemes, as well as computational complexity and efficiency of QuANNs 
over classical ANNs (Chrisley, 1995; Kak, 1995; Menneer and Narayanan, 

^This degree of freedom behaves like spin, so that the neuron’s associated qubit can 
also be approached in terms of a spin-1/2 model. 
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1995; Behrman et ai, 1996; Menneer, 1998; Ivancevic and Ivancevic, 2010; 
Gongalves, 2015). 

Mathematically, a classical ANN with a binary firing pattern can be de¬ 
fined as an artificial networked compnting system comprised of a directed 
graph with the following additional strnctnre (McCnlloch and Pitts, 1943; 
Muller et al, 1995): 

• A binary alphabet A 2 = {0,1} associated to each neuron describing the 
neural activity, with 0 corresponding to a non-firing neural state and 1 
to a firing neural state, so that the firing patterns of a neural network 
with N neurons are expressed by the set of all binary strings of length 
N: A 2 = {S 1 S 2 ...SN : Sk e A2,k = 1,2,...,A^}; 

• A real-valued weight associated with each neural link, expressing the 
strength and type of neural connection; 

• A transfer function which determines the state transition of the neuron 
and that depends upon the state of its incident neurons, the weight 
associated with each incoming neural links and an activation threshold 
that can be specific for each neuron. 

A quantum version of ANNs, on the other hand, can be defined as a directed 
graph with a networked quantum computing structure, such that (Gongalves, 
2015): 

• To each neuron is associated a two-dimensional Hilbert Space 7^2 spanned 
by the computational basis B 2 = {|0), |1)}, where |0), |1) are ket vec¬ 
tors (in Dirac’s bra-ket notation for Quantum Mechanics’ vector-based 
formalism using Hilbert space^, where |0) encodes a non-firing neural 
dynamics and |1) encodes a firing neural dynamics; 

• To a neural network, comprised of N neurons, is associated the tensor 
product of N copies of 7^2, so that the neural network’s Hilbert space 

^We use the vector representation convention introduced by Dirac (1967) for Hilbert 
spaces, assumed and used extensively in Quantum Mechanics. In this case, a ket vector, 
represented as |a), is a column vector of complex numbers while a bra vector, represented 
as (a|, is the conjugate transpose of |a), that is: (a| = |a)h The Hilbert space inner 
product is represented as (|a) , |6)) = {a\b). The outer product is, in turn, given by |a) {b\. 

A projection operator corresponds to an operator of the form Pa = |a) (a| which acts on 
any ket \b) as Pa \b) = {a\b) |o). 
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is the space spanned by the basis = ||s) : s G which 
encodes all the alternative firing patterns of the neurons; 

• The general neural configuration state of the neural network is charac¬ 
terized by a normalized ket vector |'0) G expanded in the neural 
firing pattterns’ basis as: 

1^) = V'(s) |s) (9) 

sSA^ 

with the normalization condition; 

5^lV'(s)P = l (10) 

sSA^ 

• The neural network has an associated neural links state transition op¬ 
erator Ljvet such that, given an input neural state It/'m), the operator 
transforms the input state for the neural network in an output state 
It/'oni), reflecting, in this operation, the neural links for the neural net¬ 
work, so that each neuron has an associated structure of unitary oper¬ 
ators that is conditional on its input neurons; 

IV^ont) LneI (H) 

The output state of a QuANN shows, in general, complex quantum corre¬ 
lations so that the quantum dynamics of a single neuron may depend in a 
complex way on the entire neural network’s configuration (Gongalves, 2015). 
Considering the neurons rii, ...,71^ for a A^-neuron neural network, the 
operator can be expressed as a product of each neuron’s neural links operator 
following the ordered sequence ui, ...,nAr, where neuron ni is the first to be 
updated and un the last (that is, following the activation sequenc^; 

LNet = Ln---L2Li ( 12 ) 

each neuron’s neural links operator is a quantum version of an activation 
function, with the following structure for the k-th neuron: 

Lk= Y I®) ® Lk{sin) C) |s') (s'l (13) 

seA^-hs'eA^-'' 

•^For some QuANNs it is possible to consider the action of the operators conjointly and 
to introduce, in one single neural links operator, a transformation of multiple neurons’ 
states, taking advantage of parallel quantum computation (Gongalves, 2015). 
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where Sj„ is a substring, taken from the binary word ss', that matches in 
ss' the activation pattern for the input neurons of rik, under the neural net¬ 
work’s architecture, in the same order and binary sequence as it appears in 
ss', Lk{sin) is a neural links function that maps the input substring to a uni¬ 
tary operator on the two-dimensional Hilbert space 7^2, this means that, for 
different configurations of the neural network, the neural links operator for 
the k-th neuron assigns a corresponding unitary operator that depends 
upon the activation pattern of the input neurons. 

The neural links operators incorporate the local structure of neural con¬ 
nections so that there is a unitary state transition for the neuron (a quantum 
computation) conditional upon the firing pattern of its input neurons. 

Now, an arbitrary unitary operator on a single-Hilbert space 7^2 
is a member of the unitary group U(2) and can be derived from a specific 
Hamiltonian operator structure (Greiner and Muller, 2001), so that we have, 
for a QuANN, a conditional unitary state transition: 

Tfc(sin) := (14) 

where the neuron’s associated Hamiltonian operator Hsi^ is conditional on 
the input neurons’ firing pattern and given by the general structure: 

(15) 

i=i 

where h is the reduced Planck constant]^ 6{sin), oj{sin) are measured in ra¬ 
dians per second and depend upon the neural configuration for the input 
neurons, 1 is the unit operator on 7^2, the Mj(sj„) terms are the components 
of a real unit vector u(sj„) and aj are Pauli’s operator^ 

ff, = |0){l| + |l){0|=(') J) 

0-2 = -!|0) {1| +i|l) (0| = ( ^ q’ ) 

ff3 = |0)(0|-|l)(l|= g ^ 

H.054571800(13) x 

^The terms (h,/2)dj, in the Hamiltonian, are equivalent to the spin operators for a 
spin-1/2 system (Leggett, 2002). 


( 16 ) 

( 17 ) 

( 18 ) 
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Replacing Eq.(15) in Eq.(14) and expanding we obtain; 


= e 2 


cos 






1 — i sin 


Lki^in) 6 ^ 

3 

i=i 


6'(sj„)At 


( 19 ) 


where 1 = |0) (0| + |1) (1| is the unit operator on 'H 2 - The operator in Eq.(19) 
is comprised of the product of a phase transformation exp (ici;(sj„) Af/2) and a 
rotation operator defined as (Greiner and Muller, 2001; Nielsen and Chuang, 
2003): 


R\l{Sin) [^(Sin), At] 



( 20 ) 

An arbitrary single-unitary operator (a quantum logic gate on a qubit) 
can, thus, be expressed by the product (Nielsen and Chuang, 2003): 




At 


= exp 


.tj(sj^) Af 


-Ru(..„) [9(Sm)Ai] 


( 21 ) 


This means that the transfer function of classical ANNs is replaced, for 
QuANNs, by phase transforms and rotations of the neuron’s quantum state 
conditional upon the firing pattern of the input neuron^ 

Now, given an operator O on the neural network’s Hilbert space 
expanded as: 

6=5^ (22) 

SjS'sA^ 

taking the inner product between a normalized ket vector l"^) and the trans¬ 
formed vector 0\'ip) yields ; 


,0 


= 


o 


V’) = 


SjS'SA^ 


= c>s,s'V’(s')V’(s)* 

SjS'sA^ 


(23) 


®This leads to quantum correlations that reflect the neural network’s structure 
(Gongalves, 2015). 
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For Hermitian operators obeying the relation: 

O,,,' {s'|i/>> ('/'|s> i5b,s' 


( 24 ) 


given that the state vector is normalized, if this relation is verihed, then 
Eq.(23) yields a classical expectation in which the amplitudes in square mod¬ 
ulus |'0(s) p are equivalent to decision weights associated with each alternative 
value on the diagonal of the operator’s matrix representation: 



so that, for a neural network in the state \'ip), the neural activity can be 
described by the value Os,s with an associated weight of |'^(s)p. 

In the case of Econophysics, as well as Game Theory applications, one 
usually assumes that the social system tends to the alternatives in proportion 
to the corresponding decision weights, such that one can associate a proba¬ 
bility measure for the system to follow each alternative as numerically coin¬ 
cident to the corresponding decision weight. This is akin to Game Theory’s 
notion of mixed strategy, in the sense that each player can be characterized 
by a fixed mixed strategy and play probabilistically according to the mixed 
strategy’s weights. 

While the probability of a player’s behavior is zero or one after play, the 
decision weights remain the same, in the case of Game Theory this means 
that the Nash equilibrium does not change, being available as a cognitive 
strategic scheme for further plays (Nash, 1951). In applications of QuANNs 
to social systems this means that one needs to work with either an Everettian 
interpretation of Quantum Theory, or with a Bohmian interpretatioiiQ 

The Bohmian interpretation is often assumed by researchers dealing with 
Econophysics (Ghoustova, 2007a,b; Khrennikov, 2010; Haven and Khren- 
nikov, 2013), in particular, when one wishes to address the amplitudes in 
square modulus |'^(s) p in terms of economic forces linked to emergent degrees 
of freedom that tend to make the system follow certain paths probabilistically 
(a quantum-based probabilistic version of Haken’s slaving principle applied 
to economic and hnancial systems (Haken, 1977)). 

The Everettian line of interpretations has, since its initial proposal by Ev¬ 
erett (1957, 1973), been directly linked to a Gybernetics’ paradigmatic basis 

^Since these are the two lines of interpretation that do not assume a state vector 
collapse. 
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incorporating both Automata Theory and Information Theory (Gongalves, 
2015), a point that comes directly from Everett’s original work on Quantum 
Mechanics, that is further deepened by Deutsch’s work on Quantum Compu¬ 
tation (Deutsch, 1985), and, later, on Quantum Decision Theory (Deutsch, 
1999; Wallace, 2002, 2007). 

There are actually different perspectives from different authors on Ev¬ 
erett’s original proposal (Bruce, 2004). Formally, the proposal is close to 
Bohm’s, including the importance attributed to computation and to Infor¬ 
mation Theory, however, systemically, Bohm and Everett are very distinct 
in the hypotheses they raise: for Bohm the state vector is assumed to rep¬ 
resent a statistical average of an underlying information hold’s sub-quantum 
dynamics (Bohm, 1984; Bohm and Hiley, 1993), Everett (1957, 1973) as¬ 
sumes the geometry of the Hilbert space as the correct description of the 
fundamental dynamics of helds and systems. 

Considering QuANNs, under Everett’s approach, and working from the 
state vector, we can introduce the set of projection operators onto the ba¬ 
sis V = |Ps = |s) (s| : s G I where each operator has the matrix 

representation Ps,s' = Ss,s', these operators form a complete set of orthogo¬ 
nal projectors, since their sum equals the unit operator on the Hilbert space 
= 1®'^, and they are mutually exclusive, that is, the product 

of two of these operators obeys the relation Sg^s'PsPs'- 

A projection operator can represent a projective computation, by the 
neural network, of an alternative neural bring pattern for the network. The 
general state vector in Eq.(9) can, thus, be expressed as a sum of projections, 
that is, the neural network’s quantum state has a projective expression over 
each alternative neural configuration simultaneously, corresponding to a si¬ 
multaneous systemic projective activity over all alternatives: 

|V> = E A (26) 

sSA^ 

Each alternative neural configuration corresponds to an orthogonal dimension 
of the 2^ dimensional Hilbert space a dimension that is spanned by a 

corresponding basis vector in Bf^, which means that the quantum system 
(in our case, the QuANN) projects simultaneously over each (orthogonal) 
dimension of systemic activity (corresponding, in our case, to each alternative 
neural pattern) weighing each dimension. The wheight of the projection over 
a given dimension (a given pattern of systemic activity) in the system’s state 
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can be worked from a notion of norm. Using the Hilbert space’s inner product 
structure, we can work with the squared norm of the projected vector, which 
leads to: 

|A |^)f = (A 1^), A1^)) = |AA| = |^(s)p (27) 

Systemically, this last equation can be interpreted as expressing that the 
weight of the projection Pg, in the system’s projective dynamics, is equal to 
|'^(s)p. In this sense, each orthogonal dimension corresponds to a distinct 
pattern of activity that is projectively computed by the system. 

On the other hand, for a large ensemble of QuANNs with the same struc¬ 
ture and in the same state, the statistical weight associated to the projection 
operator A, expressed by the ensemble average ^A^, coincides with the pro¬ 
jection weight associated to the neural state projection A It/’), thus, 

the statistical interpretation comes directly from the projective structure of 
the state vector. Indeed, let us consider a statistical ensemble of M QuANNs 
such that each QuANN has the same number of neurons N and the same 
architecture, let us, further, assume that each neural network is characterized 
by some quantum neural state I'^fc), with k = 1, 2,..., M, the ensemble state 
can be represented by a statistical density operator; 

1 ^ 

(28) 

^ k=l 

the statistical average of an operator O on the Hilbert space is given 
by (Bransden and Joachain, 2000); 

= Tr{dp) = 

1 ^ 

s.s'eA^^ 

. M 1 ^ 

k=i k=i ^ 

for a projector on the neural basis we get the ensemble average 

^ k=l k=l 


(29) 


( 30 ) 
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Now, if all the members of the ensemble are in the same nenral state \'ipk) = 

l"^) for each /c = 1, ...,M the whole statistical weight that is placed on the 

projection coincides exactly with |'0(s)P so that the ensemble average of 

the projection coincides nnmerically with the degree to which the system 

projects over the dimension corresponding to the nenral pattern |s) (the 

2 


projetion norm), that is, there is a nnmerical coincidence between 



and ( Ps 


p 




1 

M 


M 

E 

k=l 


S = 


|l/l(s)|^ 


(31) 


Thns, an ensemble of QnANNs with the same strnctnre, characterized 
by the same qnantnm state |V’), has a statistical weight for each projection 
coincident with the norm of the projection, so that this norm has a statistical 
expression once we consider an ensemble of systems with the same strnctnre 
and characterized by the same state. This is similar to the argnment that 
is made aronnd repeated independent]^ and identically prepared experiments 
leading to a statistical distribntion that shows the markers of the nnderlying 
qnantnm dynamics, in that case, we also see a statistical ensemble marker 
(considering an ensemble of experiments with the same state vector) that 
recovers the projection norm strnctnre in the statistical distribntion. 

The experiments, in the case of hnman systems, have led to the hnding 
of the same compntational properties and projective dynamics present in 
the qnantnm systems (Bnsemeyer and Franco, 2010; Bnsemeyer and Brnza, 
2012; Wang and Bnsemeyer, 2013; Bnsemeyer and Wang, 2014), a hnding 
that comes from the statistical distribntion of the experiments. 

In an Econophysics setting, the projective dynamics can be addressed 
as a cognitive projection snch that the projection norm corresponds to the 
decision weight placed on that alternativ^ The QnANN state transition 

®In the case of QuANNs this pressuposes the non-interaction between the ensemble 
elements, appealing to a description of a statistical random sample. 

®In the quantum computational setting, under the Everettian line, the projective struc¬ 
ture for QuANNs can be considered as a computational projection such that each Hilbert 
space dimension, corresponding to a different neural pattern, is computed simultaneously 
with an associated weight (given by the norm of the projection), having a computational 
expression in the system’s quantum processing and a statistical correspondence in the 
neural activity pattern of an ensemble of QuANNs with the same structure and in the 
same state (assuming non-interaction between different ensemble elements). 
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has an implication in the projection weights, in the sense that given the 
state transition; 


It/’ont) ^Net It/’in) 


the output amplitudes are given by: 


^ ^ '4^out(^^') I®) ) 

sSA^ 




out 


s = 


E 


s'eA^ 


^Net 


(sVm) 


s'eA^ 


(32) 


(33) 


with L 7 Vet(s,s') = Lpfet s y. Eq.(33) means that the following change in 
the projections’ norms takes place; 


-)■ 


Psl'ipout) =|^o«t(s)|^ = 


Psl'ipin) =|^m(s)p-t 


^ ^ pNetip) S ')'lpin{p ) 

^'SA^ 


(34) 


The sum within the square modulus is a source of quantum interference at 
the projection norm level. 

An iterative scheme with the repeated application of the neural network 
operator leads to a sequence of quantum neural states \'ip{t)). Expanding 
the complex numbers associated to the quantum amplitudes: 

(s|^/>(t)) = (35) 


we can express the dynamical variables A{s,t) and B{s,t) in terms of a 
dynamical nonlinear state transition rule: 


A{s,t) 


n 2 


Re 




LNet{s,s')'il;{s',t - At) 




(36) 


In the case of physical systems, the projective dynamics, interpreted computationally, 
leads to a physical expression of the system at multiple dimensions of systemic activity, 
a point which was interpreted by DeWitt (1970) under the notion of Many Worlds of a 
same Universe, where each World corresponds to an entire configuration of the Universe 
matching a corresponding orthogonal dimension of an appropriate Hilbert space where 
observers and systems are correlated (entanglement). 

In the case of applications to human decision-making, the orthogonal dimensions can be 
assumed to correspond to alternative decision scenarios evaluated by the decision-maker 
and supporting his/her choice. 
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2 

(37) 

which leads to a 2^^^ system of nonlinear equations, from where it follows 
that the probability associated to a given neural bring conbguration, worked 
from the expected projection (in accordance with the ensemble average), is 
given by the sum of the two dynamical variables: 

Prob[s,t] = A{s,t) + B{s,t) (38) 

This establishes a bridge between Nonlinear Dynamical Systems Theory 
and quantum processing by QuANNs, with implications for bnancial mod¬ 
eling. Indeed, while, traditionally, in bnancial econometrics one can see the 
distinction between a stochastic process (be it linear or nonlinear) and a 
deterministic nonlinear dynamical system, in the case of QuANNs applied 
to bnancial modeling they synthesize both approaches (stochastic and de¬ 
terministic nonlinear), since the quantum state transition equations have a 
corresponding expression in a nonlinear deterministic dynamical system for 
probability measures assigned to the QuANN’s statistical description via the 
correspondence between the projection norm dynamics and the statistical 
expectation associated to the projection operator. 

The QuANNs application to bnancial modeling, thus, allows us to address 
the problem of simulating the resulting system dynamics that comes from a 
human cognition where interference patterns are found in the probabilistic 
description of human behavior. 

2.3 A Quantum Market Model 

Considering the bnancial case, a quantum regime switching model for the N 
volatility components plus the market polarization component, introduced 
in subsection 2.1, can be addressed through a Quantum Neural Automaton 
(QNA), debned as a one dimensional lattice with a QuANN associated to 
each lattice site, in this case we assume the lattice to have iV -|- 1 sites and to 
each site k, for k = 1,2,N + 1, is associated a QuANN with the following 
structure: 

ni{k) -)■ n2{k) 

\ (39) 

nsik) 


B{s,t) = 


Im LNet{s,s')'il;{s',t - At) 


'6A^ 
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The corresponding Hilbert space for each such neural network TiNetik) is 
that is l-LNet{k) = for k = 1, 2,iV + 1, with the general basis vector 
|>siS 2 <S 3 ), such that si characterizes the activity pattern of the hrst neuron 
{ni{k)), S 2 characterizes the second neuron {n 2 {k)) and S 3 characterizes the 
activity pattern of the third neuron {n^{k)). In what follows, the neuron 
n^{k) encodes the market state for the corresponding component, ni{k) en¬ 
codes the new market conditions supporting the corresponding component’s 
dynamics and n 2 {k) addresses the computation of the synchronization pat¬ 
tern between n^{k) (the market state for the component) and ni{k) (the new 
market conditions). The QNA Hilbert space T-Lqna = ®k=i'ki-Net{k) is the 
tensor product of iV-|-1 copies of the Hilbert space Assuming this struc¬ 
ture for the QNA, we now begin by addressing the local neural dynamics and 
its financial interpretation. 

2.3.1 Local Neural Dynamics 

Since the third neuron firing patterns encode the market state of the corre¬ 
sponding component, for the N volatility components, we have the neural 
network market volatility operator on 


Ov kiS 20 ) = no kiS 20 ) 

(40) 

Ov IS 1 S 2 I) = vi IS 1 S 2 I) 

(41) 


while for the market polarization component we have the neural network 
market polarization operator; 


Op S1S2O) = -1 S1S2O) 

(42) 

Op S1S2I) = +1 S1S2I) 

(43) 


Since, as dehned previously, no ^ for a volatility neural network, when the 
third neuron hres we have a high volatility state, and when it does not hre 
we have a low volatility state. For the market polarization neural network, 
when the third neuron fires we have a bullish market state and when it does 
not fire we have a bearish market state. 

Eqs.(40) to (43) show that both operators depend only on the third neu¬ 
ron’s firing pattern, which means that, using Dirac’s bra-ket notation, they 
can be expanded, respectively, as: 
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(44) 


Ov = Vo 

Op = - 


1-S1-S2O) (S1S201 j + 

,S 2 £A 2 / \si,S 2 Sj 

IS1S2O) (S1S201 j + f 

!1,S2SA2 / \si,S2S-' 


IS1S21) (S1S21I 


I-S1S21) (S1S21I 

,S2SA2 J 


] 


(45) 


Now, the neural network follows a closed loop starting at the market state 
neuron {no{k)) and ending at the market state neuron. The hnal state tran¬ 
sition amplitudes and the underlying hnancial dynamics will depend upon 
the intermediate transformations which may change the prohle of the corre¬ 
sponding component’s state transition structure. 

To address the neural dynamics and its relation with the hnancial market 
dynamics we need to introduce the neural links operators and follow the loop, 
starting at nolk) and ending at no{k). Considering, then, the hrst neural link 
nsik) —t ni{k), we introduce the following neural network operator for the 
neuron ni{k)\ 


Li = e 0 |Q^ ^Q| ^ g ^ 1 ^/^ 

- 5 GA 2 


using Eq.(19) we need to dehne the angles 6*(0), 6^(1), a;(0), a;(l) and the 
vectors u(0), u(l), we set, in this case: 


6*(0)At 71 9{l)At 

2 ~ 2 ’ 2 ~ 

a;(0)At a;(l)At tt 

2 “ 2 ~ 2 

u(0)=u(l) = (1,0,0) 
leading to the following operator structure: 






sin (0) i cos (0) 
i cos (0) sin (0) 


= i cos (0) 1 -|- sin (0) di = 
the action of the operator Li on the basis states is given by: 


i cos (0) sin (0) 
sin (0) i cos (0) 


(47) 

(48) 

(49) 

(50) 

(51) 
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Li IS 1 S 2 O) = sin (0) IS 1 S 2 O) + i cos (0) |1 — S 1 S 2 O) (52) 

Li IS 1 S 2 I) = i cos (0) IS 1 S 2 I) + sin (0) |1 - S 1 S 2 I) (53) 

The operator Li can be considered in terms of a quantum regime switch¬ 
ing model, such that if the market state neuron ns{k) is not bring, then, 
sin (0) is the amplitude associated to the alternative where the neuron ni{k) 
does not change state, while i cos (0) is the amplitude associated to the alter¬ 
native where the neuron ni{k) changes state, on the other hand, if the neuron 
ns{k) is firing the role of the amplitudes flip; icos{(j)) is associated with the 
alternative where the neuron ni{k) does not change state and sin (0) is the 
amplitude associated with the alternative where the neuron ni{k) changes 
state. 

Before considering the financial implications of this dynamics, it is nec¬ 
essary to address the rest of the network, because the final dynamics and its 
financial implications can only be fully addressed at the end of the cycle. As 
we will see, the end result will be a quantum computation-based selection 
process of adaptive rules regarding market expectations and the processing 
of how financial news may support trading decisions affecting market polar¬ 
ization and market volume. 

Proceeding, then, with the neural links, the second neuron to be activated 
is n 2 {k), which receives an input from the two neurons ni{k) and n^{k), 
this neuron will play a key role in the selection of adaptive rules regarding 
the relation between trading prohles and financial news, a point that we 
will return to when the final neural network state transition is analyzed. 
Following the quantum circuit framework, the second neuron is transformed 
conditionally on the states of the two neurons ni{k) and n^{k), in accordance 
with the neural links ni{k) —)■ n 2 {k) rislk), the corresponding neural links 
operator is given by: 

1/2 = ^ l-s) (s| ® ® |s') (s'l (54) 


When the input neurons have synchronized firing patterns, the rotation and 
phase transformation angles are set to; 


0(OO)At _ e{ll)At 
2 ~ 2 ~ 

a;(00)Af _ u{ll)At 
2 ~ 2 ~ 


(55) 

(56) 
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which means that the operators reduce to: 




( 57 ) 


that is, the second neuron remains in the same state when the input neu¬ 
rons {ni{k) and n^{k)) exhibit a synchronized bring pattern (no rotation nor 
phase transformation takes place). When the input neurons do not exhibit a 
synchronized bring pattern, the rotation and phase transformation is set by 
the following parameters: 


0(Ol)At e{ 10 )At n 

2 “ 2 “2 

(58) 

a;(01)At a;(10)At tt 

2 “ 2 “2 

(59) 

u(01) = u(10) = (1,0,0) 

(60) 

which leads to: 


g-|Athoi _ 

(61) 

thus, the action of L2 on each basis state is such that: 


L2 SS2S) = SS2S) 

(62) 

L2 SS2I - s) = si - S2I - s) 

(63) 


that is, the neuron n 2 {k) does not change state when the two neurons ni{k) 
and n^lk) have the same bring pattern, and bips state when the two neu¬ 
rons have dibering bring patterns (this is equivalent to a controlled negation 
quantum circuit). 

Now, to close the cycle, and before addressing the bnal dynamics and 
its bnancial interpretation, we have to address, brst, the third link n 2 {k) —)■ 
n^{k). In this case, we also introduce a controlled-negation circuit, so that 
the corresponding operator is: 

I /3 = ^ k) (s| ® |0) (0| ® ^ |s') (s'l ® |1) (1| ® (64) 

seA2 s'eA2 


6'(0)At 6*(l)At TT 

2 ~ ’ 2 ~ 2 

a;(0)Af a;(l)At tt 

2 “ ’ 2 “ 2 


(65) 

( 66 ) 
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u(l) = (1,0,0) 

(67) 

leading to: 


e-^Atho ^ ^ 

( 68 ) 

SO that the basis states transform as: 


L3 kiOsa) = IsiOss) 

(69) 

L 3 \S 1 IS 3 ) = |sill - S 3 ) 

(70) 


these equations show that the neuron rislk) changes state when the second 
neuron is bring and does not change state when the second neuron is not br¬ 
ing. The neural network operator is the product of the three operators, 
that is; 


Lnci — L3L2L1 


(71) 


The following table shows the results of the action of the neural network 
operator on each basis state. 


Basis States 

LNet s) = L 3 L 2 L 1 s) 

|000) 

sin (0) 000) -I- i cos (0) 111) 

0 

0 

i cos (0) 010) 4- sin (0) 101) 

0 

0 

sin (0) oil) -I- i cos (0) 100) 

1011) 

i cos (0) 001) -I- sin (0) 110) 

1100) 

i cos (0) 000) -I- sin (0) 111) 

1101) 

sin (0) 010) -1- i cos (0) 101) 

1110) 

i cos (0) oil) -I- sin (0) 100) 

1111) 

sin (0) 001) -1- i cos (0) 110) 


Table 1: Neural network operator’s action on the basis states. 

From a bnancial perspective, table 1 synthesizes two adaptive rules, one 
in which the new market state for the component follows the new market con¬ 
ditions underlying the corresponding component’s dynamics (neurons’ ni{k) 
and ns{k) show a neural reinforcement dynamics), and another in which the 
new market state is contrarian with respect to the new market conditions un¬ 
derlying the corresponding component’s dynamics (neurons’ ni{k) and nslk) 
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show a neural inhibitory dynamics). These are two basic rules regarding ex¬ 
pectation formation from new data: the decision to follow the new data or 
not. 

In the hrst case, and taking as example a volatility component, the market 
is driven by an expectation of continuance of market conditions, so that, for 
instance, if market conditions are favorable to a high volatility state (neuron 
ni{k) is bring), then, the new market state follows the market conditions and 
n^{k) fires, corresponding to high volatility. 

On the other hand, still under the hrst adaptive rule, if market conditions 
are unfavorable to a high volatility state (neuron ni{k) is not firing), then, 
the new market state follows the market conditions and n^{k) does not hre, 
corresponding to low volatility. 

The resulting adaptive rule corresponds, thus, to a follow the news rule. 
Likewise, if we consider, instead, the market polarization component, the 
follow the news rule means that if the new market conditions support a 
bullish market sentiment, then, the market becomes bullish and if the new 
market conditions support a bearish market sentiment, then, the market 
becomes bearish. 

The second adaptive rule is the reverse, expectations are that the new 
market conditions will not hold, and the market does the opposite from the 
news, expecting speculative gains. 

The hrst adaptive rule is implemented when the second neuron is not 
bring, while the second rule is implemented when the second neuron is br¬ 
ing. Thus, the bring of the second neuron is a dynamical component that 
simulates a market change in its expectation and trading proble, so that, 
for the neural conbgurations {|000), |001), |100), |101)}, the state transition 
for the market component’s dynamics is driven by the hrst adaptive rule, 
while, for the neural conbgurations {|010), |011), |110), the market 

component’s dynamics is driven by the second adaptive rule. 

While neuron n 2 (/c)’s bring pattern determines the selection of a follow the 
news rule, the combination of bring patterns of the three neurons determines 
the quantum amplitudes for the market state transitions. Thus, when the 
neuron n 2 {k) is not bring, if the initial market conditions are aligned with 
the initial market state, then; sin (0) is the amplitude associated with the 
alternative in which ni{k) and n^lk) transition to a not bring state and 
icos (0) is the amplitude associated with the alternative in which ni{k) and 
ns{k) transition to a bring state. 

For a market volatility component, this means that a transition to a high 
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volatility state, supported by market conditions, has an associated quantum 
amplitude of icos(0), while a market transition to a low volatility, state 
supported by market conditions, has an associated amplitude of sin (0). The 
role of these amplitudes switches when ni{k) and n^{k) are not initially 
aligned. 

When the neuron n 2 {k) is bring, the transition amplitudes to bring/non- 
bring states follow the same pattern as above for neuron ni{k) but reverse 
the pattern for neuron n‘i{k) because the new market conditions’ neuron and 
the market state neuron transition to a non-aligned state (the market is con¬ 
trarian with respect to the news), so that, if ni{k) and n^{k) are initially 
aligned, sin (0) is the amplitude associated with a transition to the state 
where ni{k) is not bring and n^{k) is bring, while, if ni{k) and n^{k) are not 
initially aligned, the amplitude associated with such a transition is icos (0). 
The roles of the amplitudes, thus, depend upon the way in which the market 
adapts to new information and the previous conbguration of market condi¬ 
tions and market state. 

As expected, the market conditions and the market state neurons are al¬ 
ways entangled, which means that, in each case, the market state ebectively 
becomes like a measurement apparatus of the market conditions, the entan¬ 
glement proble can, however, be aligned {follow the news rule, based on an 
expectation of sustainability of the new market conditions) or non-aligned 
{contrarian rule, based on the expectation of reversal of the new market 
conditions). 

Thus, in the model, the quantum neural dynamics models a market that 
processes the information on the market conditions implementing a stan¬ 
dard quantum measurement, but the proble of that quantum measurement 
depends upon the expectations regarding the news (leading to diberent en¬ 
tanglement probles). 

The bnal dynamics for the market component results from the iterative 
application of the operator L^et for each trading round, leading to state tran¬ 
sition between the adaptive rules and, thus, between the market states. Con¬ 
sidering a sequence of neural states for the market component’s associated 
neural network |0(A;,f)), the state transition resulting from the dynamical 
rule is given by: 


\'ip{k,t)) = L]^et\'ip{k,t - At)) (72) 

leads to the following update rule for the quantum amplitudes ipk (as per the 
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general Eq.(33)): 

= X] (s |Ljvet| s'^ (s'|^/>(A;,t - At)) = 

(73) 

= ^ L^etis, s')^/Jk{s', t - At) 

using Table I’s results, in conjunction with this last eqution, we obtain the 
following transition table for the quantum amplitudes: 


New Amplitudes 

■^^(000, t) = sin (0) 0fc(OOO, t — At) + i cos (0) 0^(100, t — At) 
0fc(OOl, t) = i cos (0) 0fc(Oll, t — At) + sin (0) 0^(111, t — At) 
0fc(OlO, t) = i cos (0) 0fc(OOl, t — At) + sin (0) 0^(101, t — At) 
0fc(Oll, t) = sin (0) 0fc(OlO, t — At) + i cos (0) 0^(110, t — At) 
0fc(lOO, t) = i cos (0) 0A:(O1O, t — At) + sin (0) 0^(110, t — At) 
0fc(lOl, t) = sin (0) 0fc(OOl, t — At) + i cos (0) 0^(101, t — At) 
0fc(llO, t) = sin (0) 0fc(Oll, t — At) + i cos (0) 0^(111, t — At) 
0fc(lll, t) = i cos (0) 0fc(OOO, t — At) + sin (0) 0^(100, t — At) 


Table 2: Update of the quantum amplitudes for a single market component. 

Taking into account this general neural dynamics for each component 
we can now piece it all together to address the market state and resulting 
hnancial dynamics. 

2.3.2 Financial Market Dynamics 

To address the full market dynamics we need to recover the QNA. For each 
trading round, the quantum state associated with the market dynamics is 
given by the QNA state dehned as the tensor product of the lattice site’s 
neural networks’ states, that is, by the tensor product of each component’s 
neural network state: 

Af+l 

I^W) = 0 ^)) = S2, •••, Sv+1, t) |Sl, S2,..., Sjv+l) (74) 

k=l 

where the quantum amplitudes 'h(si,S 2 , ...,S 7 v+i,t) are given by: 

7V+1 

T(Si,S 2 , ...,SAr+i,t) = Jj0fc(Sfc,t) (75) 

fc=l 
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with ipk being the amplitudes associated with the lattice site /c’s neural net¬ 
work. 

For the N volatility components we can introduce a corresponding volatil¬ 
ity operator on the QNA Hilbert space: 

Ok = ® ® ( 76 ) 

with k = where, as before, 1®™ denotes m-tensor product of the 

unit operator on ^2 and Oy is the volatility operator dehned in Eqs.(40) 
and (41). Similarly, for the market polarization operator, we write: 

On+i = ® Op (77) 

where Op is the market polarization operator dehned in Eqs.(42) and (43). 
In this way, the returns’ dynamical variable dehned in Eq.(8) is replaced, 
in the Quantum Econophysics setting, by a quantum operator on the QNA 
Hilbert space dehned as: 

, N+l 

R=TYlOk (78) 

^ k=l 

For each basis state of the QNA Hilbert space, the returns operator has an 
eigenvalue given by the corresponding hnancial market returns: 

R |si, S2, ..., Sat+i) = i?(si, S2, ..., Sat+i) |Si, S2, ..., Sat+i) (79) 

with the eigenvalues R(si, S2, ..., Sat+i) given by: 

1 ^ 

R(si, S2, ..., Sat+i) = ^ t^fe(Sfc) ■ aAr+i(sAT+i) (80) 

^ k=l 

where Vk(sk) = Vq if the binary string s*, G A 2 ends in 0 (nz(k) is not bring) 
and nfc(sfc) = vi if the binary string Sk ends in 1 (n^(k) is firing), similarly 
(TAT+i(sAr+i) = —1 if Sat+1 G A| ends in 0 (ns(N + 1) is not firing) and 
o'N+ii^N+i) = 1 if Sat+i ends in 1 (ns(N + 1) is bring). 

The dynamical rule that comes from the neural networks’ quantum com¬ 
putation leads to the market state transition for each trading round: 

7V+1 

l-ipit)) = (^ LNet\'ip(k,t - At)) (81) 

k=l 
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leading to the expected value for the returns; 

/?(Si, S 2 , Stv+i) |d'(Si, S 2 , Sjv+l, t)|^ (82) 

si,S2,.--,Siv+i 

SO that the market tends to the alternative -R(si, S 2 ,s^r+i) with an associ¬ 
ated probability of |\['(si, S 2 ,s^v+i, t)f. 

The following hgure shows a market simulation on Python 3.4. In the 
simulations, the initial state for each component is taken from a randomly 
chosen U(2) gate applied to each neuron with uniform probability over U(2) . 
The hgure shows the markers of hnancial turbulence in the returns, including 
volatility bursts and jumps. 



Figure 1; Simulation of the hnancial returns for sin^0 = 0.6, Vq = 0.7, 
A = 1000, 20 components (19 volatility components plus 1 polarization com¬ 
ponent). The hgure shows 2000 data points of a 2100 data points simulation 
with the hrst 100 points removed for transients. 

The main parameters that determine the market prohle with regards to 
turbulence is vo and the number of components, the turbulence prohle does 
not change much with respect to the rotation angle 0. Indeed, as shown in 
the following table, the estimated kurtosi^for diherent simulations with 20 
components tends to decrease as Vq rises. For vq = 0.9 we no longer hnd 
excess kurtosis, the turbulence markers being lost. This approach to low 
turbulence is progressive as uq is raised from 0.8 to 0.9, such that that the 
price jumps tend to become less severe and less frequent, and the volatility 

^°The Fisher kurtosis is used in the statistical analysis of the model’s outputs. 
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bursts tend to disappear, as shown in Figure 2, in which vq = 0.9 with the 
rest of the parameters used in Figure I’s simulation being left unchanged. 



sin^ (p = 0.4 

sin^ (p = 0.5 

sin^ p = 0.6 

o 

II 

o 

585.4546 

1336.6726 

1923.3159 

Vo = 0.5 

778.0387 

1876.3810 

783.0852 

O 

II 

o 

1015.5296 

473.4505 

383.9775 

o 

II 

o 

77.6054 

49.5857 

56.8335 

00 

o 

II 

o 

6.8827 

20.7037 

5.6217 

no = 0.9 

-0.8538 

-1.2335 

-0.9277 


Table 3; Kurtosis values for different values of sin^ 0 and uq- The other 
parameters are: A = 1000, 20 components (19 volatility components plus 
1 polarization component), the Kurtosis coefficient was calculated on 5000 
sample data points of a 5100 data points simulation with the hrst 100 data 
points removed for transients. 



Figure 2; Simulation of hnancial returns for sin^ (f) = 0.6, Vq = 0.9, A = 1000, 
20 components (19 volatility components plus 1 polarization component). 
The hgure shows 2000 data points of a 2100 data points simulation with the 
hrst 100 points removed for transients. 

The model, thus, captures different market prohles; as the parameter uq 
increases from 0.8 to 0.9 the simulations tend to approach a lower tail risk 
dynamics, with a greater approximation to the classical Gaussian returns’ 
probability model ocurring for uq near 0.87, the following table shows this 
approximation with the kurtosis values and Jarque-Bera test for normality, 
as the value of vq is increased. 
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Kurtosis 

JB Statistic 

p-value 

VO = 0.85 

1.2911 

353.7757 

0.0 

uo = 0.86 

0.3596 

46.5669 

7.7289e-ll 

uo = 0.87 

0.1746 

6.3179 

0.0425 

uo = 0.88 

-0.0160 

76.5797 

0.0 

uo = 0.89 

-0.6790 

106.5589 

0.0 

Vo = 0.9 

-0.8223 

143.5983 

0.0 


Table 4: Kurtosis values and Jarque-Bera test of normality for different values 
Uo- The other parameters are: sin^ 0 = 0.6, A = 1000, 20 components (19 
volatility components plus 1 polarization component), the Kurtosis coefficient 
was calculated on 5000 sample data points of a 5100 data points simulation 
with the hrst 100 data points removed for transients. 

As shown in table 4 for every value of Vq the Jarque-Bera’s null hypoth¬ 
esis is rejected at 1% signihcance level except for uq = 0.87. It is important 
to stress however, that although simulated returns distribution can approx¬ 
imate the Gaussian distribution, this approximation is not robust, different 
simulations for the same parameters may show deviations from the Gaussian 
distribution. 

The following table shows examples of simulations for different values of 
the rotation angle 0, with uq = 0.87, the null hypothesis of Jarque-Bera’s 
test is not reject, at a 1% signihcance, for sin^ 0 = 0.1, 0.3, 0.6, with sin^ 0 = 
0.1, 0.3 as the only cases in which it is not rejected for 5% signihcance, 
and sin^0 = 0.3 as the only case in which it is not rejected also for a 10% 
signihcance. 
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sin^ 0 

Kurtosis 

JB Statistic 

p-value 

0.1 

0.1285 

4.8697 

0.0876 

0.2 

0.6385 

84.4538 

0.0 

0.3 

0.0947 

3.8794 

0.1437 

0.4 

0.3169 

35.1954 

2.2773e-08 

0.5 

0.9213 

565.5135 

0.0 

0.6 

0.1746 

6.3179 

0.0425 

0.7 

0.0841 

29.6076 

3.7221e-07 

0.8 

0.4975 

65.859 

4.9960e-15 

0.9 

0.1775 

31.3417 

1.5640e-07 


Table 5: Kurtosis values for different values (j). The other parameters are: 
Uo = 0.87 A = 1000, 20 components (19 volatility components plus 1 polar¬ 
ization component), the Kurtosis coefficient was calculated on 5000 sample 
data points of a 5100 data points simulation with the hrst 100 data points 
removed for transients. 


These results may, however, depend, as stated previously, upon the simu¬ 
lation, other simulations may show the null hypothesis being rejected for the 
same parameters, which means that the Gaussian distribution depends upon 
the sample path and is not a dynamically hxed probability law that can be 
assumed to hold indehnitely. 

The general tail risk pattern, on the other hand, is more robust than 
the Gaussian approximation, in the sense that as uq approaches 0.9 and 
for vq > 0.9, the market loses the turbulence prohle with the jumps and 
volatility changes becoming less frequent and the kurtosis becoming less and 
less leptokurtic, leading to lower tail risk, the market returns eventually 
fluctuate randomly around a narrow band. 

Underlying the complex behavior of the simulated market returns, is the 
probability dynamics that comes from the neural network’s iterative scheme 
shown in table 2. Gonsidering Eqs.(35) to (37) and combining with table 
2’s results we get, in this case, sixteen nonlinear dynamical equations of the 
general form; 


Ak{s,t) 


^/Ak{s',t - At) sin (0) - ^/Bk{s",t - At) cos (0) 


(83) 




— At) sin (0) -I- ^/Ak{s",t — At) cos (0) 


(84) 
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with s, s', s" G A 2 and s' 7 ^ s", so that the probability dynamics that come 
from the neural network’s evolution can be addressed by a nonlinear map 
with sixteen dynamical variables satisfying the normalization rule: 

^ Afc(s,t) + ^5fc(s',t) = 1 (85) 

S 

with the probability of the neural conhguration s being given by the sum; 

Probk[s, t] = Ak{s, t) + i?fc(s, t) (86) 


so that the probability distribution for the neural conhgurations is a function 
of a sixteen dimensional nonlinear map on a hypersphere of unit radius (due 
to the normalization condition). 

If we expand the squares in Eqs.(83) and (84) we get: 


Afc(s, t) = Afc(s', t — At) sin^ (0) + Bk{s”, t — At) cos^ (0) — 
— \/Ak{s',t — At)Bk{s”,t — At) sin (20) 
Bk{s, t) = i?fc(s', t — At) sin^ ( 0 ) + Afc(s", t — At) cos^ (0) + 
+ A/fifc(sVf^^^t))4^0s0T^^^t) sin ( 20 ) 
which leads to the following expansion for the probability; 


(87) 

( 88 ) 


Probk[s,t] = 
= Probk[s', t — At] sin^ (0) + 
+Probk[s", t - At] cos^ (0) + 
+ a/ Bk{s', t - At)Ak{s", t - At) sin (20) - 
— ^/Ak{s\ t — At)Bk{s”, t — At) sin (20) 


the quantum interference terms (that correspond to the square root terms 
multiplied by sin (20) in Eq.(89)) have an expression, at the probability level, 
that can be approached in terms of a classical nonlinear dynamical system 
for the probabilities. 

In the classical nonlinear dynamics representation, each hnancial returns 
component’s stochastic dynamics has a probability measure that updates at 
each trading round with a deterministic nonlinear update rule, this estab¬ 
lishes the bridge between the stochastic process and the nonlinear determin¬ 
istic dynamical systems modeling of hnancial dynamics: the neural network’s 
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quantum dynamics leads to a nonlinear deterministic dynamics in the prob¬ 
abilities. 

A question that may be raised regards the transition from the determin¬ 
istic nonlinear map to a noisy nonlinear map, from the hnancial perspective 
this makes sense since external stochastic factors may affect the financial 
system. A possible solution for this might be to allow the rotation angle 0 
to change, so that instead of a hxed value of 0 we replace it by a random 
variable in Eqs.(83) and (84) so that we get a stochastic nonlinear dy¬ 
namical system. The introduction of a random (j)k{t) implies that we are no 
longer dealing with a hxed unitary operator structure for the QuANN but, 
instead, work with a quantum neural state transition with a random compo¬ 
nent in the Hamiltonian, that is, the unitary gates of Eqs.(50) and (51) are 
now stochastic unitary gates: 


fAiitop) _ / 

^ sin(0fc(t)) 

icos(0fc(t)) \ 

(90) 

“ 1 

y icos(0fc(t)) 

sin(0fc(t)) J 

fAiifip) _ / 

^ i cos {4>k{t)) 

sin((/>k(t)) \ 

(91) 

“ \ 

y sin(0fc(t)) 

icos((/>k(t)) ) 


Thus, a stochastic nonlinear map is induced by the quantum noisy gates in 
the QuANN’s state transition rule, coming from a stochastic Hamiltonian. 
The following hgure shows the simulation results for 

/ 1 

dkit) = arcsin , 

^ ^ ^ VVl + e-2/3^Ht) 

with Zkit) ~ iV(0,1), which leads to; 

sin' (0fc(t)) = ^ ^ (93) 

„-2l3zk (t) 

cos^Mt))= (94) 

the logistic function present in Eqs.(93) and (94) is also widely used in clas¬ 
sical ANNs for the activation probability and leaves room for expansion of 
connections to Statistical Mechanics (Muller and Strickland, 1995). If we 
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replace in Eq.(89) we get the nonlinear stochastic equations for the proba¬ 
bilities: 

Probk[s,t] = 

Probk[s', t — At] 

X _|_ g-2/32fe(t) 

Probk[s'\t — 

1 -f e-2/3zfc(i) (95) 

1 4- g-2/3zfc(q 
g (b 
X g-2/32:fe(b 


0015 
0010 
0005 
OOOO 
- 0.005 
- 0.010 
- 0.015 
- 0.020 

Figure 3; Simulation of hnancial returns for noisy gates with (5 = 2.0, 
uo = 0.9, A = 1000, 80 components (79 volatility components plus 1 po¬ 
larization component). The hgure shows 4000 data points of a 4100 data 
points simulation with the hrst 100 points removed for transients. 

The hgure 3 shows the occurrence of price jumps and clustering volatility, 
the turbulence in this case is linked to the high number of components (rather 
than to the noisy gates). Indeed, as the following two tables show, the 
noisy gates do not have a strong effect on the transition from leptokurtic to 
platikurtic distributions, both for low and high values of /I, it is the number 
of components that have a stronger impact on market prohle, as seen in 
table 7 for the case of Uq = 0.88 which for the simulation with (3 = 2 was 
close enough to the Gaussian distribution for the non-rejection of the null 
hypothesis of the Jarque-Bera test at a 10% signihcance level. 
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(3 = 0.01 

Kurtosis 

JB p-value 

II 

to 

Kurtosis 

JB p-value 

uo = 0.86 

0.9171 

0.0 

Vo = 0.86 

0.7097 

0.0 

Vq = 0.87 

0.0345 

1.0518e-07 

Vo = 0.87 

0.1786 

0.0018 

Vo = 0.88 

-0.3647 

2.8422e-13 

Vo = 0.88 

-0.0856 

0.3119 

Vq = 0.89 

-0.7142 

0.0 

Vo = 0.89 

-0.6205 

0.0 

o 

II 

o 

-0.6810 

0.0 

O 

II 

o 

-0.8840 

0.0 


Table 6 ; Kurtosis values and Jarque-Bera test of normality p-values for dif¬ 
ferent simulations with variying vq and noisy unitary gates. The other pa¬ 
rameters are: jS = 0.01 (left table) and /9 = 2 (right table), A = 1000, 20 
components (19 volatility components plus 1 polarization component), the 
kurtosis coefficient was calculated on 5000 sample data points of a 5100 data 
points simulation with the hrst 100 data points removed for transients. 


d = 0.01 

Kurtosis 

p-value 

3 = 2 

Kurtosis 

p-value 

N + l = 10 

-1.3408 

0.0 

N + l = 10 

-1.3752 

0.0 

N + l = 20 

-0.3647 

2.8422e-13 

N + 1 = 20 

-0.0856 

0.3119 

N +1 = 30 

1.6227 

0.0 

N + l = 30 

2.2008 

0.0 

N + 1 = 40 

2.5130 

0.0 

N + 1 = 40 

3.4597 

0.0 

N + l = 50 

10.1352 

0.0 

N + 1 = 50 

15.9802 

0.0 


Table 7: Kurtosis values and Jarque-Bera test of normality for different values 
of the number of components (iV -|- 1) and noisy unitary gates. The other 
parameters are: Uq = 0.88 = 0.01 (left table) and (3 = 2 (right table), 

A = 1000, the kurtosis coefficient was calculated on 5000 sample data points 
of a 5100 data points simulation with the hrst 100 data points removed for 
transients. 

Indeed, the number of components shows a strong effect, as can be seen 
in hgure 3, which uses Vq = 0.9 and in table 7, that shows the transition from 
platikurtic to leptokurtic for large values of the components, for both a low 
and a high value of (3. 
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3 Finance, Nonlinear Stochastic Dynamics and 
Quantum Artificial Intelligence 

From the early onset of development of Econophysics, some form of non¬ 
linear stochastic dynamics has been considered to be present in hnancial 
market dynamics. A major example being Vaga’s work that addressed ex¬ 
plicitly different probability distributions corresponding to different (classi¬ 
cal) Hamiltonian conditions (Vaga, 1990). The major point that markets 
make transitions between different regimes and different probability distri¬ 
butions was key to Vaga’s market theory. On the other hand, the multifractal 
multiplicative cascades (Mandelbrot et ah, 1997) introduced multiplicative 
stochastic processes as sources of market turbulence. 

While a division line is drawn in regards to nonlinear deterministic pro¬ 
cesses versus nonlinear stochastic processes, the possible combination of both 
might provide an intermediate approach, combining adaptive market dynam¬ 
ics and stochastic factors affecting market behavior. 

As the previous section model shows, when recurrent QuANNs are applied 
to hnancial modeling, the nonlinear deterministic dynamics and the non¬ 
linear stochastic processes result directly from the quantum computational 
structure, in the sense that; while the iterative computation of a QuANN 
results from the linear conditional unitary state transition, the corresponding 
probabilities, due to the square modulus rule for addressing the probabilities 
associated to different neural bring patterns, leads to a nonlinear update rule 
for the probabilities themselves, which means that the market behavior will 
show an interference effect at the probability level expressable in terms of a 
classical nonlinear map, thus, while the system follows a stochastic dynamics, 
the probabilities are updated nonlinearly. 

This is a direct consequence of Quantum Cognitive Science that comes 
from human decision analysis, which shows that the nonlinear update in 
probabilities, leading to non-additive decision weights may be computation¬ 
ally approached from linear unitary quantum computation on an appropriate 
Hilbert space. Stochastic factors in the nonlinear update of probabilities can 
also be introduced through unitary noise in the neural network’s computation 
through stochastic Hamiltonians. 

Although QuAI and QuANN theory are still on their early stages, they 
provide a bridge between major lines of research on hnancial dynamics and 
risk modeling including; nonlinear deterministic and stochastic dynamics ap- 
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plied to financial modeling, Cognitive Science and computational foundations 
of Financial Theory. Future research on QuANNs dynamics may thus serve 
as a relevant tool to link different approaches that characterized the different 
lines of research on Econophysics-based Finance. 
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